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Abstract. In the heat-bath Glauber dynamics for the Ising model on 

the lattice, physicists believe that the spectral gap of the continuous- 
time chain exhibits the following behavior. For some critical inverse- 
temperature /3c, the inverse-gap is 0(1) for (3 < /3c, polynomial in the 
surface area for = f3c and exponential in it for 13 > 0c. This has been 
proved for I? except at criticality. So far, the only underlying geometry 
where the critical behavior has been confirmed is the complete graph. 
Recently, the dynamics for the Ising model on a regular tree, also known 
as the Bethe lattice, has been intensively studied. The facts that the 
inverse-gap is bounded for /3 < /3c and exponential for /3 > /3c were 
established, where 0c is the critical spin-glass parameter, and the tree- 
height h plays the role of the surface area. 

In this work, we complete the picture for the inverse-gap of the Ising 
model on the 6-ary tree, by showing that it is indeed polynomial in h 
at criticality. The degree of our polynomial bound does not depend 
on 6, and furthermore, this result holds under any boundary condition. 
We also obtain analogous bounds for the mixing-time of the chain. In 
addition, we study the near critical behavior, and show that for /3 > /3c, 
the inverse-gap and mixing-time are both exp[©((/3 — /3c)/i)]. 



The Ising Model on a finite grapli G = (V, E) with inverse-temperature 
/? > and no external magnetic field is defined as follows. Its set of possible 

configurations is O = {±1}^, where each cj G assigns positive or negative 
spins to the vertices of G. In the free boundary case, the probability that 
the system is at a given configuration a is given by the Gibbs distribution 



where Z(l3) is the partition function. In the presence of a boundary condition 
T G {±1}^^ (that fixes the spins of some subset dV C F of the sites), we let 
IJ'q{o') denote the Gibbs measure conditioned on a agreeing with r on dV. 

The heat-bath Glauber dynamics for the Ising model on G is the Markov 
chain with the following transition rule: At each step, a vertex is chosen 
uniformly at random, and its spin is updated according to /Iq conditioned 
on the spins of all the other vertices. It is easy to verify that this chain is 
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reversible with respect to the Gibbs distribution fiQ. The continuous-time 
version of the dynamics associates each site with an independent Poisson 
clock of unit rate, determining the update times of this site as above (note 
that the continuous dynamics is |y| times faster than the discrete dynamics). 

The spectral-gap of a reversible discrete-time chain, denoted by gap, is 
1 — A, where A is the largest nontrivial eigenvalue of its transition kernel. 
The spectral-gap of the continuous-time process is defined analogously via 
the spectrum of its generator, and in the special case of Glauber dynamics 
for |Ug, this gap is precisely \V\ times the discrete-time gap. This parameter 
governs the rate of convergence to equilibrium in L^(//g). 

In the classical Ising model, the underlying geometry is the d-dimensional 
lattice, and there is a critical inverse-temperature /3c where the static Gibbs 
measure exhibits a phase transition with respect to long-range correlations 
between spins. While the main focus of the physics community is on critical 
behavior (see the 20 volumes of [14]), so far, most of the rigorous mathe- 
matical analysis was confined to the non-critical regimes. 

Supported by many experiments and studies in the theory of dynamical 
critical phenomena, physicists believe that the spectral-gap of the continuous- 
time dynamics on lattices has the following critical slowing down behavior 
(e.g., [16,20,25,38]): At high temperatures (/? < (3c) the inverse-gap is 0(1), 
at the critical (3c it is polynomial in the surface-area and at low temperatures 
it is exponential in it. This is known for except at the critical (3c, and 
establishing the order of the gap at criticality seems extremely challenging. 
In fact, the only underlying geometry, where the critical behavior of the 
spectral-gap has been fully established, is the complete graph (see [9]). 

The important case of the Ising model on a regular tree, known as the 
Bethe lattice, has been intensively studied (e.g., [3-7,15,17,18,23,27,31]). 
On this canonical example of a non- amenable graph (one whose boundary is 
proportional to its volume), the model exhibits a rich behavior. For example, 
it has two distinct critical inverse-temperatures: one for uniqueness of the 
Gibbs state, and another for the purity of the free-boundary state. The 
latter, Pc, coincides with the critical spin-glass parameter. 

As we later describe, previous results on the Ising model on a regular tree 
imply that the correct parameter to play the role of the surface-area is the 
tree-height h: It was shown that the inverse-gap is 0(1) for (3 < (3c and 
exponential in h for (3 > (3c, yet its critical behavior remained unknown. 

In this work, we complete the picture of the spectral-gap of the dynamics 
for the critical Ising model on a regular tree, by establishing that it is indeed 
polynomial in h. Furthermore, this holds under any boundary condition, 
and an analogous result is obtained for the L} (total-variation) mixing time, 
denoted by tMix (formally defined in Subsection 2.1). 



MIXING OF CRITICAL ISING MODEL ON TREES 



3 



Theorem 1. Fix b > 2 and let (3c = arctanh(l/-v/6) denote the critical 
inverse-temperature for the Ising model on the b-ary tree of height h. Then 
there exists some constant c > independent ofb, so that the following holds: 
For any boundary condition t, the continuous-time Glauber dynamics for the 
above critical Ising model satisfies gap~"^ < tMix = 0{h^). 

One of the main obstacles in proving the above result is the arbitrary 
boundary condition, due to which the spin system loses its symmetry (and 
the task of analyzing the dynamics becomes considerably more involved). 
Note that, although boundary conditions are believed to only accelerate the 
mixing of the dynamics, even tracking the effect of the (symmetric) all-plus 
boundary on lattices for /? > /3c is a formidable open problem (see [26]). 

In light of the above theorem and the known fact that the inverse-gap is 
exponential in h at low temperatures {(3 > (5c fixed), it is natural to ask how 
the transition between these two phases occurs, and in particular, what the 
critical exponent oi (5 — (3c is. This is answered by the following theorem, 
which establishes that log(gap~^) x (/? — (3c)h + log/i for small (3 — (3c- 
Moreover, this result also holds for P = (3c -\- o{l), and thus pinpoints the 
transition to a polynomial inverse-gap at (3 — (3c'^ 

Theorem 2. For some £q > 0, any b > 2 fixed and all (3c < (3 < (3c + Eq, 
where (3c = arctanh(l/\/6) is the critical spin-glass parameter, the following 
holds: The continuous-time Glauber dynamics for the Ising model on a b-ary 
tree with inverse-temperature (3 and free boundary satisfies 

gap^^ = exp [Q {{(3 — (3c)h)] otherwise. 

Furthermore, both upper bounds hold under any boundary condition t, and 
(1.1) remains valid if ga.p~^ is replaced by tuix- 

In the above theorem and in what follows, the notation / = Q{g) stands 
for / = 0(5) and g = 0{f). 

Finally, our results include new lower bounds on the critical inverse-gap 
and the total- variation mixing-time (see Theorem 3). The lower bound 
on gap^^ refutes a conjecture of [3], according to which the continuous-time 
inverse-gap is linear in h. Our lower bound on ^mix is of independent interest: 
Although in our setting the ratio between ^mix and gap~^ is at most poly- 
logarithmic in n, the number of sites, we were able to provide a lower bound 
of order logn on this ratio without resorting to eigenfunction analysis. 

1.1. Background. The thoroughly studied question of whether the free 
boundary state is pure (or extremal) in the Ising model on the Bethe lattice 
can be formulated as follows: Does the effect that a typical boundary has 
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on the spin at the root vanish as the size of the tree tends to infinity? It is 
well-known that one can sample a configuration for the tree according to the 
Gibbs distribution with free boundary by propagating spins along the tree 
(from a site to its children) with an appropriate bias (see Subsection 2.2 for 
details). Hence, the above question is equivalent to asking wether the spin 
at the root can be reconstructed from its leaves, and as such has applications 
in Information Theory and Phylogeny (see [15] for further details). 

In sufficiently high temperatures, there is a unique Gibbs state for the 
Ising model on a 6-ary tree {b > 2), hence in particular the free boundary 
state is pure. The phase-transition with respect to the uniqueness of the 
Gibbs distribution occurs at the inverse-temperature /3„ = arctanh(l/6), as 
established in 1974 by Preston [34]. 

In [6], the authors studied the critical spin-glass model on the Bethe lattice 
(see also [5, 7]), i.e., the Ising model with a boundary of i.i.d. uniform spins. 
Following that work, it was finally shown in [4] that the phase-transition in 
the free-boundary extremality has the same critical inverse-temperature as 
in the spin-glass model, /3c = arctanh(l/\/5). That is, the free-boundary 
state is pure iff /? < (3c- This was later reproved in [17, 18]. 

The inverse-gap of the Glauber dynamics for the Ising model on a graph 
G was related in [3] to the cut-width of the graph, ^(G), defined as follows: 
It is the minimum integer m, such that for some labeling of the vertices 
{vi, . . . , Vn} and any k £ [n], there are at most m edges between {vi, . . . , v^} 
and {vk+i, ■ ■ ■ ,Vn}- The authors of [3] proved that for any bounded degree 
graph G, the continuous-time gap satisfies gap~^ = exp[0(^(G)/?)]. 

Recalling the aforementioned picture of the phase-transition of the gap, 
this supports the claim that the cut-width is the correct extension of the 
surface- area to general graphs. One can easily verify that for (the d- 
dimensional box of side-length L) the cut-width has the same order as the 
surface-area L'^~^, while for a regular tree of height h it is of order h. 

Indeed, for the Ising model on a 6-ary tree with h levels and free boundary, 
it was shown in [3] that the inverse-gap is 0(1) for all (3 < /3c, whereas for 
(3 > (3c '^^ satisfies loggap^^ x h (with constants that depend on b and (3). 
The behavior of the gap at criticality was left as an open problem: it is 
proved in [3] that the critical gap~^ is at least linear in h and conjectured 
that this is tight. A weaker conjecture of [3] states that gap~^ = exp(o(/i)). 

Further results on the dynamics were obtained in [27], showing that the 
log-Sobolev constant as (defined in Section 2) is uniformly bounded away 
from zero for /3 < /3c in the free-boundary case, as well as for any P under the 
all-plus boundary condition. While this implies that gap~^ = 0(1) in these 
regimes, it sheds no new light on the behavior of the parameters gap, as in 
our setting of the critical Ising model on a regular tree with free-boundary. 
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1.2. The critical inverse-gap and mixing-time. Theorems 1,2, stated 
above, establish that on a regular tree of height h, the critical and near- 
critical continuous-time gap~^ and tuix are polynomial in h. In particular, 
this confirms the conjecture of [3] that the critical inverse-gap is exp(o(/i)). 

Moreover, our upper bounds hold for any boundary condition, while 
matching the behavior of the free-boundary case: Indeed, in this case the 
critical inverse-gap is polynomial in h (as [3] showed it is at least linear), 
and for j3 — 13c > small we do have that log(gap^^) x (/? — (3c)h. For 
comparison, recall that under the all-plus boundary condition, [27] showed 
that gap~^ = 0(1) at all temperatures. 

We next address the conjecture of [3] that the critical inverse-gap is in 
fact linear in h. The proof that the critical gap~^ has order at least h uses 
the same argument that gives a tight lower bound at high temperatures: 
Applying the Dirichlet form (see Subsection 2.4) to the sum of spins at 
the leaves as a test-function. Hence, the idea behind the above conjecture 
is that the sum of spins at the boundary (that can be thought of as the 
magnetization) approximates the second eigenfunction also for /3 = (3c- 

The following theorem refutes this conjecture, and en-route also implies 
that Os = o(l) at criticality. In addition, this theorem provides a nontrivial 
lower bound on tuix that separates it from gap~^ (thus far, our bounds in 
Theorems 1,2 applied to both parameters as one). 

Theorem 3. Fix b > 2 and let (3c = arctanh(l/\/6) be the critical inverse- 
temperature for the Ising model on the b-ary tree with n vertices. Then the 
corresponding discrete-time Glauber dynamics with free boundary satisfies: 

gap"^ > cin(logn)^ , (1.2) 
W > C2n{\ognf , (1.3) 

for some ci, C2 > 0. Furthermore, tuix > cgap~^ logn for some c > 0. 

Indeed, the above theorem implies that in continuous-time, gap^^ has 
order at least h'^ and tMix has order at- least h"^, where h is again the height 
of the tree. By known facts on the log-Sobolev constant (see Section 2, 
Corollary 2.4), in our setting we have ^mix = 0{a~^h), and it then follows 
that Os = 0{h-^) = o(l). 

We note that by related results on the log-Sobolev constant, it follows 
that in the Ising model on a regular tree, for any temperature and with 
any boundary condition we have ^mix = 0(gap~^ log^ n). In light of this, 
establishing a lower bound of order log n on the ratio between tMix and gap~^ 
is quite delicate (e.g., proving such a bound usually involves constructing a 
distinguishing statistic via a suitable eigenfunction (Wilson's method [39])). 
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1.3. Techniques and proof ideas. To prove the main theorem, our gen- 
eral approach is a recursive analysis of the spectral-gap via an appropriate 
block-dynamics (roughly put, multiple sites comprising a block are updated 
simultaneously in each step of this dynamics; see Subsection 2.5 for a formal 
definition) . This provides an estimate of the spectral-gap of the single-site 
dynamics in terms of those of the individual blocks and the block-dynamics 
chain itself (see [25]). However, as opposed to most applications of the 
block-dynamics method, where the blocks are of relatively small size, in our 
setting we must partition a tree of height h to subtrees of height linear in h. 
This imposes arbitrary boundary conditions on the individual blocks, and 
highly complicates the analysis of the block-dynamics chain. 

In order to estimate the gap of the block-dynamics chain, we apply the 
method of Decomposition of Markov chains, introduced in [19] (for details 
on this method see Subsection 2.6). Combining this method with a few 
other ideas (such as establishing contraction and controlling the external 
field in certain chains), the proof of Theorem 1 is reduced into the following 
spatial-mixing/reconstruction type problem. Consider the procedure, where 
we assign the spins of the boundary given the value at the root of the tree, 
then reconstruct the root from the values at the boundary. The key quantity 
required in our setting is the difference in the expected outcome of the root, 
comparing the cases where its initial spin was either positive or negative. 

This quantity was studied by [31] in the free-boundary case, where it was 
related to capacity- type parameters of the tree (see [15] for a related result 
corresponding to the high temperature regime). Unfortunately, in our case 
we have an arbitrary boundary condition, imposed by the block-dynamics. 
This eliminates the symmetry of the system, which was a crucial part of the 
arguments of [31]. The most delicate step in the proof of Theorem 1 is the 
extension of these results of [31] to any boundary condition. This is achieved 
by carefully tracking down the effect of the boundary on the expected re- 
construction result in each site, combined with correlation inequalities and 
an analytical study of the corresponding log-likelihood-ratio function. 

The lower bound on the critical inverse-gap reflects the change in the 
structure of the dominant eigenfunctions between high and low tempera- 
tures. At high temperatures, the sum of spins on the boundary gives the 
correct order of the gap. At low temperatures, a useful lower bound on 
gap^-*^ was shown in [3] via the recursive-majority function (intuitively, this 
reflects the behavior at the root: Although this spin may occasionally flip its 
value, at low temperature it constantly tends to revert to its biased state). 
Our results show that at criticality, a lower bound improving upon that of [3] 
is obtained by essentially merging the above two functions into a weighted 
sum of spins, where the weight of a spin is determined by its tree level. 
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To establish a lower bound on tuix of order gap^ /i, we consider a cer- 
tain speed-up version of the dynamics: a block-dynamics, whose blocks are 
a mixture of singletons and large subtrees. The key ingredient here is the 
Censoring Inequality of Peres and Winkler [33], that shows that this dy- 
namics indeed mixes as fast as the usual (single-site) one. We then consider 
a series of modified versions of this dynamics, and study their mixing with 
respect to the total-variation and Hellinger distances. In the end, we arrive 
at a product chain, whose components are each the single-site dynamics on 
a subtree of height linear in h. This latter chain provides the required lower 
bound on ^Mix* 



1.4. Organization. The rest of this paper is organized as follows. Section 2 
contains several preliminary facts and definitions. In Section 3 we prove a 
spatial-mixing type result on the critical and near-critical Ising model on a 
tree with an arbitrary boundary condition. This then serves as one of the 
key ingredients in the proof of the main result. Theorem 1, which appears 
in Section 4. In Section 5 we prove Theorem 3, providing the lower bounds 
for the critical inverse-gap and mixing-time. Section 6 contains the proof 
of Theorem 2, addressing the near-critical behavior of gap~^ and tMix- The 
final section. Section 7, is devoted to concluding remarks and open problems. 



2. Preliminaries 

2.1. Total- variation mixing. Let (Xt) be an aperiodic irreducible Markov 
chain on a finite state space Q, with stationary distribution tt. For any two 
distributions (p, ip on 0,, the total-variation distance of (p and tp is defined as 




The (worst-case) total-variation mixing-time of {Xt), denoted by tMix(e) for 
< e < 1, is defined to be 

W(e) = min \ t : max ||Px(-'^'t G •) " ^IItv < ^ \ , 

where denotes the probability given that Xq = x. As it is easy and 
well known (cf., e.g., [1, Chapter 4]) that the spectral-gap of (Xt) satisfies 
gap~^ < tuix (l/e)i it will be convenient to use the abbreviation 

^MIX = ^MIX (1/g) • 

Analogously, for a continuous-time chain on 0, with heat-kernel Ht, we define 
tMix as the minimum t such that max^.gn \\Ht{x, •) — vtHtv < 
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2.2. The Ising model on trees. When the underlying geometry of the 
Ising model is a tree with free boundary condition, the Gibbs measure has 
a natural constructive representation. This appears in the following well 
known claim (see, e.g., [15] for more details). 

Claim 2.1. Consider the Ising model on a tree T rooted at p with free 
boundary condition and at the inverse-temperature (3. For all e G E{T), 
let ?7e S {il} be i.i.d. random variables with P(r/e = 1) = (1 + tanh/3)/2. 
Furthermore, let a{p) be a uniform spin, independent of {r]e}, and for v ^ p, 

a{v) = a{p) 7]e ) where V{p,v) is the simple path from p to v. 

Then the distribution of the resulting a is the corresponding Gibbs measure. 

In light of the above claim, one is able to sample a configuration accord- 
ing to Gibbs distribution on a tree with free boundary condition using the 
following simple scheme: Assign a uniform spin at the root p, then scan 
the tree from top to bottom, successively assigning each site with a spin 
according to the value at its parent. More precisely, a vertex is assigned the 
same spin as its parent with probability (1 + tanh/?)/2, and the opposite 
one otherwise. Equivalently, a vertex inherits the spin of its parent with 
probability tanh /3, and otherwise it receives an independent uniform spin. 
Finally, for the conditional Gibbs distribution given a plus spin at the root 
p, we assign p a plus spin rather than a uniform spin, and carry on as above. 

However, notice that the above does not hold for the Ising model in the 
presence of a boundary condition, which may impose a different external 
influence on different sites. 

2.3. L^-capacity. The authors of [31] studied certain spatial mixing prop- 
erties of the Ising model on trees (with free or all-plus boundary conditions), 
and related them to the L^-capacity of the underlying tree. In Section 3, 
we extend some of the results of [31] to the (highly asymmetric) case of a 
tree with an arbitrary boundary condition, and relate a certain "decay of 
correlation" property to the L^-capacity of the tree, defined as follows. 

Let T be a tree rooted at p, denote its leaves by 9T, and throughout the 
paper, write (n, v) G -E(T) for the directed edge between a vertex u and its 
child V. We further define dist(ii,u) as the length (in edges) of the simple 
path connecting u and v in T. 

For each e G E{T), assign the resistance Re > to the edge e. We say 
that a non-negative function / : E(T) — > M is a flow on T if the following 
holds for all (u, v) G E{T) with v dT: 

f{u,v) = f{v,w) , 
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that is, the incoming flow equals the outgoing flow on each internal vertex 
V in T. For any How /, define its strength \ f\ and voltage V{f) by 

1/1= E f{p,v),V{f)^supl E f{e)Re:wGdT\ 

where V{p, w) denotes the simple path from p to w. Given these definitions, 
we now define the L^-capacity cap2(T) to be 

cap2(T) = sup{|/| : / is a flow with V{f) < 1} . 

For results on the L^-capacity of general networks (and more generally, L^- 
capacities, where the expression f{e)Re in the above definition of V{f) is 
replaced by its (p— 1) power), as part of Discrete Nonlinear Potential Theory, 
cf., e.g., [28], [36], [37] and the references therein. 

For our proofs, we will use the well-known fact that the L^-capacity of 
the tree T is precisely the effective conductance between the root p and the 
leaves (?T, denoted by Cefr(p ^ dT). See, e.g., [24] for further information 
on electrical networks. 

2.4. Spectral gap and log-Sobolev constant. Our bound on the mixing 
time of Glauber dynamics for the Ising model on trees will be derived from a 
recursive analysis of the spectral gap of this chain. This analysis uses spatial- 
mixing type results (and their relation to the above mentioned capacity) 
as a building block. We next describe how the mixing-time can be bounded 
via the spectral-gap in our setting. 

The spectral gap and log-Sobolev constant of a reversible Markov chain 
with stationary distribution vr are given by the following Dirichlet forms 
(see, e.g., [1, Chapter 3,8]): 

where 

£{f) = ((/ - P)f, /), = ^ E [/(^) - /(y)]' <^)n^^ y) ' (2-2) 

Ent,(/)=E,(/2log(/VlE./')) . (2.3) 

By bounding the log-Sobolev constant, one may obtain remarkably sharp up- 
per bounds on the mixing-time: cf., e.g., [10-13,35]. The following result 
of Diaconis and Saloff-Coste [12, Theorem 3.7] (see also [35, Corollary 2.2.7]) 
demonstrates this powerful method; its next formulation for discrete-time 
appears in [1, Chapter 8]. As we are interested in total- variation mixing, we 
write this bound in terms of tMix, though it in fact holds also for the (larger) 
mixing-time. 
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Theorem 2.2 ([12], [35], reformulated). For any reversible finite Markov 
chain with stationary distribution vr, 

Wx(l/e) < \a-^ (loglog(l/7r*) +4) , 

where tt* = min^^ 7r(x) . 

We can then apply a result of [27], which provides a useful bound on Og in 
terms of gap in our setting, and obtain an upper bound on the mixing-time. 

Theorem 2.3 ([27, Theorem 5.7]). There exists some c > such that the 
Ising model on the b-ary tree with n vertices satisfies ag > c • gap/ log n. 

Note that the proof of the last theorem holds for any (5 and under any 
boundary condition. Combining Theorems 2.2 and 2.3, and noticing that 
TT* > 2~" exp(— 2/3n) (as there are 2" configurations, and the ratio be- 
tween the maximum and minimum probability of a configuration is at most 
exp(2/3n)), we obtain the following corollary: 

Corollary 2.4. The Glauber dynamics for the Ising model on a b-ary tree 
with n vertices satisfies ^mix = O (a^-'^logn) = O (gap^^log^n) for any [) 
and any boundary condition. 

The above corollary reduces the task of obtaining an upper bound for the 
mixing-time into establishing a suitable lower bound on the spectral gap. 
This will be achieved using a block dynamics analysis. 

2.5. From single site dynamics to block dynamics. Consider a cover of 
y by a collection of subsets . . . , i^fc}, which we will refer to as "blocks" . 
The block dynamics corresponding to . . . , iJ^ is the Markov chain, where 
at each step a uniformly chosen block is updated according to the stationary 
distribution given the rest of the system. That is, the entire set of spins of 
the chosen block is updated simultaneously, whereas all other spins remain 
unchanged. One can verify that the block dynamics is reversible with respect 
to the Gibbs distribution 

Recall that, given a subset of the sites U <Z V , a boundary condition r/ 
imposed on U is the restriction of the sites = V \ U to all agree with 77 
throughout the dynamics, i.e., only sites in U are considered for updates. It 
will sometimes be useful to consider ry G (rather than a configuration of 
the sites U^), in which case only its restriction to U'^ is accounted for. 

The following theorem shows the useful connection between the single- 
site dynamics and the block dynamics. This theorem appears in [25] in 
a more general setting, and following is its reformulation for the special 
case of Glauber dynamics for the Ising model on a finite graph with an 
arbitrary boundary condition. Though the original theorem is stated for 
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the continuous-time dynamics, its proof naturally extends to the discrete- 
time case; we provide its details for completeness. 

Proposition 2.5 ([25, Proposition 3.4], restated). Consider the discrete 
time Glauber dynamics on a h-ary tree with boundary condition ij. Let gap^ 
be the spectral-gap of the single-site dynamics on a subset U C V of the 
sites, and gap^ be the spectral-gap of the block dynamics corresponding to 
Bi, . . . , Bk, an arbitrary cover of a vertex set W <ZV . The following holds: 

gap^ > TT^gap^ inf inf |Bi|gapg^ (sup #{i : Bi 3 x}] 

Proof. Let P denote the transition kernel of the above Glauber dynamics. 
Defining 

g = infinf |Sj|gapg. , 
the Dirichlet form (2.1) gives that, for any function /, 

Combining this with definition (2.2) of (?(•), 

On the other hand, definition (2.2) again implies that 



= E/^^(V')^ EE E l^%icTm\P^^{a,a^ma) - fia^f 

< I sup #{i : B, B x}Y,l^W E \W\PU<^,a^)[f{a) - f^a^)]' 

= \W\ sup #{i : Bi 3 x}£^{f) , 
xew 

where is the configuration obtained from a by flipping the spin at x, and 
we used the fact that 

|Si|Pg^((T, C7^) = \ W\P^{a, a"") for any i G [k] and x e Bi . 

Altogether, we obtain that 

SUf) < sup #{i : Bi B x}S^^{f) . 
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Recalling that the single-site dynamics and the block-dynamics have the 
same stationary measure, 

Var;^(/) Var^(/) " ^^^^ 
(where we again applied inequality (2.1)), thus 

Var4.(/) \W\ V^gp^ ) 

The proof is now completed by choosing / to be the eigenfunction that 
corresponds to the second eigenvalue of (achieving gap^), with a final 
application of (2.1). ■ 

The above proposition can be applied, as part of the spectral gap analysis, 
to reduce the size of the base graph (though with an arbitrary boundary 
condition), provided that one can estimate the gap of the corresponding 
block dynamics chain. 

2.6. Decomposition of Markov chains. In order to bound the spectral 
gap of the block dynamics, we require a result of [19], which analyzes the 
spectral gap of a Markov chain via its decomposition into a 'projection chain 
and a restriction chain. 

Consider an ergodic Markov chain on a finite state space with transition 
kernel P : 17 x $7 — > [0, 1] and stationary distribution vr : J7 — > [0, 1]. We 
assume that the Markov chain is time-reversible^ that is to say, it satisfies 
the detailed balance condition 

7r(j;)P(x, y) = 7r(y)P(y, x) for all x, y S . 

Let O = r^oU. . .ur2m-i be a decomposition of the state space into m disjoint 
sets. Writing [m] = {0, . . . , m — 1}, we define vf : [m] — > [0, 1] as 

7f(i) = ■ni^li) = ^ 7r(x) 
and define P : [m] x [m] — > [0, 1] to be 

The Markov chain on the state space [m] whose transition kernel is P is called 
the projection chain, induced by the partition Oq, . . . ,r2m-i- It is easy to 
verify that, as the original Markov chain is reversible with respect to vr, the 
projection chain is reversible with respect to the stationary distribution vf. 
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In addition, each J7j induces a restriction chain, whose transition kernel 
Pj : rjj X — > [0, 1] is given by 



P{x,y), lix^y, 

Again, the restriction chain inherits its reversibihty from the original chain, 
and has a stationary measure vTj, which is simply vr restricted to O^: 

TTi{x) = 7r(a;)/7f(i) for all x G ilj . 

In most applications, the projection chain and the different restriction chains 
are all irreducible, and thus the various stationary distributions vf and 
TTo, . . . , TTm-i are all unique. 

The following result provides a lower bound on the spectral gap of the 
original Markov chain given its above described decomposition: 

Theorem 2.6 ([19, Theorem 1]). Let P he the transition kernel of a finite 
reversible Markov chain, and let gap denote its spectral gap. Consider the 
decomposition of the chain into a projection chain and m restriction chains, 
and denote their corresponding spectral gaps by gap and gapQ, . . . , gap^_]^. 
Define 

gaPmin = gapi , 7 = max max V P{x,y) . 

Then gap, the spectral gap of the original Markov chain, satisfies: 

^.r. > sap gap • gap^i^ 
gap > — A — — - — z — . 
3 37 + gap 

The main part of Section 4 will be devoted to the analysis of the projection 
chain, in an effort to bound the spectral gap of our block dynamics via the 
above theorem. 



3. Spatial mixing of Ising model on trees 

In this section, we will establish a spatial-mixing type result for the Ising 
model on a general (not necessarily regular) finite tree under an arbitrary 
boundary condition. This result (namely. Proposition 3.1) will later serve as 
the main ingredient in the proof of Theorem 1 (see Section 4). Throughout 
this section, let /3 > be an arbitrary inverse-temperature and 9 = tanh /?. 

We begin with a few notations. Let T be a tree rooted at p with a 
boundary condition r G {±1}^-^ on its leaves, and /i"^ be the corresponding 
Gibbs measure. 

For any v £ T, denote by T„ the subtree of T containing v and its all 
descendants. In addition, for any B <Z A <Z T and a G {±1}"^, denote by 
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aB the restriction of a to the sites of B. We then write /ujj for the Gibbs 
measure on the subtree given the boundary tqTv ■ 

Consider T C T \ dT, a subtree of T that contains the root p, and write 
Tv = Tv n T. Similar to the above definitions for T, we denote by fi^ the 
Gibbs measure on T given the boundary condition ^ G {±1}^"^, and let fA 
be the Gibbs measure on T„ given the boundary ^qj, . 

The following two measures are the conditional distributions of /xJJ on the 
boundary of the subtree Ty given the spin at its root v. 

Qt^O = ^;(^at„ = W„ I ^(^) = 1) for ^ e {±1}'''^ ' 

QviO = /^;(^af„ = W„ I ^i^) = -1) for ^ e i^^f^ ■ 



We can now state the main result of this section, which addresses the prob- 
lem of reconstructing the spin at the root of the tree from its boundary. 

Proposition 3.1. Let T be as above, let < d < ^ and define 

Then there exists an absolute constant k > ^ such that 

cap2(f) 

- K(l-e) ' 

where the resistances are assigned to be R(u,v) — 6~'^^^^^'^P''"\ Furthermore, 
this also holds for any external field h £M on the root p. 

To prove the above theorem, we consider the notion of the log-likelihood 
ratio at a vertex v with respect to Ty given the boundary : 



4=fog rr ^ ' K (3.1) 

\pl{a{v) = -1)/ 

as well as the following quantity, analogous to A (defined in Proposition 
3.1): 

my = J xldQy - J xldQ~ . (3.2) 

As we will later explain, > for any v £ T, and we seek an upper 
bound on this quantity. One of the main results of [31] was such an esti- 
mate for the case of free boundary condition, yet in our setting we have an 
arbitrary boundary condition (adding a considerable amount of difficulties 
to the analysis). The following theorem extends the upper bound on rup to 
any boundary; to avoid confusion, we formulate this bound in terms of the 
same absolute constant k given in Proposition 3.1. 
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Theorem 3.2. Let T and rup he as above, and let < 9 < |. Then there 
exists an absolute constant k > tm such that 




where the resistances are assigned to be R(u,v) 



0—2 dist(p,?;) 



Proof of Theorem 3.2. As mentioned above, the novelty (and also the 
main challenge) in the result stated in Theorem 3.2 is the presence of the 
arbitrary boundary condition r, which eliminates most of the symmetry that 
one has in the free boundary case. Note that this symmetry was a crucial 
ingredient in the proof of [31] for the free boundary case (namely, in that 
case and Q~ are naturally symmetric). 

In order to tackle this obstacle, we need to track down the precise influence 
of the boundary condition r on each vertex v G T. We then incorporate this 
information in the recursive analysis that appeared (in a slightly different 
form) in [23] . This enables us to relate the recursion relation of the m^-s to 
that of the L'^-capacity. 

The following quantity captures the above mentioned influence of r on a 
given vertex v € T: 



Notice that x* has a similar form to xi (defined in (3.1)), and is essentially 
the log-likelihood ratio at v induced by the boundary condition r. The 
quantity xt is then the log-likelihood ratio that in addition considers the 
extra constraints imposed by ^. Also note that a free boundary condition 
corresponds to the case where x* = for all v £ T. 

To witness the effect of x*, consider the probabilities of propagating a 
spin from a parent v to its child w, formally defined by 



we define pi ^{—1,1) and pJJ ^(— 1, — 1) analogously. The next simple lemma 
shows the relation between x* and these probabilities. 

Lemma 3.3. The following holds for any {v, w) G T: 




(3.3) 



pI,w(.'^, 1) = l^l{(^{w) = 1 I a{v) = 1) , 
-1) = l^l{cr{w) = -1 I cr{v) = 1) ; 



pUl,l)-pl^J-l,l) = Die 
where D*^ = (cosh /5) V (cosh^ /5 + cosh2(x;^/2) - l) . 
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Proof. Recalling definition (3.3) of x*, we can translate the boundary con- 
dition T into an external field on the vertex w when studying the 
distribution of its spin. Hence, 



e' 



e2/3 _ e-2/3 



3^™ + e + e^^ + e 
cosh^ /? 



tanh /? 



cosh^ P + cosh"'(a;^/2) - 1 
as required. ■ 

Remark 3.4. In the free boundary case, we have p„,«,(l, 1) —Pv,w{—^, 1) = d- 
For a boundary condition r, the coefficient < < 1 represents the 
contribution of this boundary to the propagation probability. 

We now turn our attention to my. As mentioned before, the fact that 
ruy > follows from its definition (3.2). Indeed, the monotonicity of the 
Ising model implies that the measure stochastically dominates the mea- 
sure Q~ (with respect to the natural partial order on the configurations of 
dTy). For instance, it is easy to see this by propagating 1 and —1 spins from 
the root to the bottom, and applying a monotone coupling on these two pro- 
cesses. Finally, xf, is monotone increasing in ^ (again by the monotonicity 
of the Ising model), thus m„ > 0. 

The first step in establishing the recursion relation of (that would 
lead to the desired upper bound) would be to relate to some quantities 
associated with its children, as stated next. 

Lemma 3.5. For any w € T \ dT, we have that 



niv = ^ 



where 

, , /cosh(x/2) -h6'sinh(2;/2)\ , 

Proof. We need the following well-known lemma, that appeared in [23] in a 
slightly different form; see also [2] and [31, Lemma 4.1]. 

Lemma 3.6 ([23], [31] (reformulated)). Let f be as in (3.4). For all v G 

f\df and ^ G {±1}^'^, the following holds: xl = Y.w:iv,w)ef fi^^)- 
According to this lemma, we obtain that 

= E ( / fi^i)dQtiO - J fi^i)dQ~{0) ■ (3.5) 

W.{v,w)S:T 
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Noting that Xw is actually a function of S^gf , we get that 

f{xi)dQt{0 = I /(xi)d(pj;,^(i,i)g+(0 + p;,^(i,-i)Q-(0) , 

and similarly, we have 

f{oci)dQ-io = J /(xi)d(p;^(-i,i)gi(o+p;^(-i,-i)g^(0) • 

Combining these two equalities, we deduce that 

f{xi)dQt{i)- j f{xi)dQ-{i) 

= (p;,.(l,l)-p;,.(-l,l))( j f{xi)dQUC)- J f{xi)dQ-{0) 

= eDl[j f{xi)dQl{i) - j f{xi)dQ-{i)) , (3.6) 

where in the last inequality we applied Lemma 3.3. Plugging (3.6) into (3.5) 
now completes the proof of the lemma. ■ 

Observe that in the free boundary case, QtiO — Qv{~0 for any ^. 
Unfortunately, the presence of the boundary r breaks this symmetry, causing 
the distributions and to become skewed. Nevertheless, we can still 
relate these two distributions through the help of formulated by the 

following lemma. 

Lemma 3.7. For v G T, let Qj, be the following distribution over {ztl}^"^.- 

QviO = QUO = fJ-l {^^9f, = iaf. 
Then for all ^ G {±1}'^-^, we have 

i * 

QUO - QviO = C: (tanh ^ - tanh ^)q.{0 , 
where C* = 2 cosh2(j;;/2). 

Proof. It is clear from the definitions of Q^, Q~ and that 
^+,,, Qv{Ol^U<^{y) = 1 I 1 + tanh(4/2) 



IJl{a{v) = l) 1 + tanh(a;*/2) 

_ QviOl^lHv) = -l\0 _ l-tanh(x|/2) 
~ iiZ{a{v) = -1) " 1 - tanh(x*/2) "^"^^^ ' 
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Hence, a straightforward calculation gives that 

^ _ 2(tanh(x|/2) -tanh(x;/2)) 

Qv iO - Qv (0 = (1 + tanh(x*/2))(l - tanh(x*/2)) ^'^^^^ 

= 2cosh2(^)(tanh^-tanh^)Q,(e) , 

as required. ■ 

The following technical lemma will allow us to combine Lemmas 3.5, 3.7 
and obtain an upper bound on in terms of {ruyj : {v,w) G T}. Note 
that the constant n mentioned next is in fact the absolute constant k in the 
statement of Theorem 3.2. 

Lemma 3.8. Let f be defined as in (3.4) for some < 9 < |. Then 

\f{x) - f{y)\ <2f{\x - y\/2) for any x,y , (3.7) 

and there exists a universal constant k > such that for any two constants 
Ci, C2 > 1 with C2 > 1 + (^Ci - 1) (1 - 9'^) and any 6 > we have 

f{5) {1 + 4k{1 - 9)Ci6tanh{6/2)) < 0266 . (3.8) 
Proof. We first show (3.7). Put 6 = \x — y\, and define 

hi6)=snp\fit + S)-fit)\ . 
t 

We claim that 

h{6) = f{5/2) - f{-5/2) = 2/(5/2) . (3.9) 

The second equality follows from the fact that f{x) is an odd function. To 
establish the first equality, a straightforward calculation gives that 

fix) 



l + (l-^2)sinh2(x/2) ' 

and it follows that f'{x) is an even non-negative function which is decreasing 
in X > 0. The following simple claim therefore immediately implies (3.9): 

Claim 3.9. Let g{t) > be an even function that is decreasing in t > 0. 
Then G{t) = fl; g{x)dx has G{t + 5)- G{t) < 2G[5/2) for any t and 6 > 0. 

Proof. Fix 6 > and define F{t) as follows: 

F{t) = G{t + 6)- G{t) . 

We therefore have F'{t) = g{t + 5) — g{t). Noticing that 

\t + 5\>\t\ ift>-| 
|t-|-(5| < \t\ otherwise ' 
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the assumption on g now gives that F'{t) < while i > — | and otherwise 
F'{t) > 0. Altogether, we deduce that 



Fit) < F{ 



G( 



G{ 



2G( 



2) ' 



as required. 



It remains to show that (3.8) holds for some k = k{6q) > 0. Clearly, it 
suffices to establish this statement for C2 = [l + (5C1 ~ l) (1 ~ ^^)] V I and 
any Ci > 1. Rearranging (3.8), we are interested in a lower bound on 



inf 

e<eo , t>o , Ci>i 



i + (f -i)(i-02)\ vij et-f{t) 



4Ci{l - e)tf{t) tanh(|) 

First, consider the case 1 < Ci < 2. We then have C2 
expression being minimized in (3.10) takes the form: 



(3.10) 
1, and the 



9t-f{t) 



4Ci(l -0)t/(i)tanh(^ 



> 



et-fit) 



[l-9)tf{t)tanh{i 



T{t,9) , 



where the inequality is by our assumption that Ci < 2. We therefore have 
that infe<gp ^ t>o r(t, 9) minimizes (3.10) for Ci < 2, and will next show that 
this is also the case for Ci > 2 under a certain condition. Indeed, letting 



9it,9,Ci) 



[l + (f 



1 a 



')] ot - fit) 



4Ci(l - 6')t/(t)tanh(t/2) 
it is easy to verify that the following holds: 



dg 



m 



dCi 4^2(1 -0)t/(t)tanh(§) ' 

hence g is increasing in Ci for every 9, t such that f{t) > 9^t. Therefore, 

g{t, 9, Ci) > g{t, 9, 2) = r(t, 9) for ah t, 9 such that f{t) > 9h . 

Before analyzing T{t, 9), we will treat the values of 9, t such that f{t) < 9'^t. 
Assume that the case, and notice that the numerator of g then satisfies 



1 + 



iyi-e^)]9t-fit) 



> 



1 + 



2 



1 (1 



9t 



and thereby the dependency on Ci vanishes: 



g{t,9,Ci) > 



9{l + 9) 



4(1 - 9)tf{t) tanh(t/2) 8f{t) tanh(t/2) 
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Since both tanh(i/2) and f{t) are monotone increasing in t and are bounded 
from above by 1 and log ( j respectively, we get 



where the second inequality is by the fact that log(l + x) < x for any x > 0, 
and the last inequality follows by the assumption 9 < ^ . 

It thus remains to establish a uniform lower bound on T(t,9). In what 
follows, our choice of constants was a compromise between simplicity and 
the quality of the lower bound, and we note that one can easily choose 
constants that are slightly more optimal. 

Assume first that ^ > ^o ^ for some to be defined later. Notice that 

and so f{t, 9) is strictly increasing in 9 for any t > 0. Since 

9t-f{t) 



> 



8(1 -^)t/(t)tanh(|) 

Ot-fit) t-f{t) 



8(1 - 9o)tf{t) tanh(|) 8(1 - 9o)tf{t) tanh(i) 

we have that T is monotone decreasing in 9 for any such t, and therefore 
8(i3eg)f (t), where f is defined as follows: 

t(t) = with respect to ^ = | . (3.12) 

^' t/(t,0)tanh(|) ^ ^ ^ ' 

Recall that the Taylor expansion of /(i, 9) around is - -^^^t-" + 0{r) . 
It is easy to verify that for ^ = | this function satisfies 

f{t, 9)<9t- for = I and any < t < 3 . 

Adding the fact that tanh(x) < x for all x > 0, we immediately obtain that 

f > = ^> 7:7^ for all < t < 3 . 

- 20t{9t){t/2) 10 20 

On the other hand, for t > 3 we can use the uniform upper bounds of 1 and 
log(Y3|) for tanh(t/2) and f{t) respectively, and gain that 

- 0t-log(i±f) 9 1 1 

r > ^r^f^ = TTfl > 7^ for alH > 3 . 

- tlog(i±f) log(f±f) t-20 

Altogether, as F > 8(j^i^f , we can conclude that T > [160(1 — ^o)]~^- 



MIXING OF CRITICAL ISING MODEL ON TREES 



21 



Note that the trivial choice of 60 = already provides a uniform lower 
bound of for T (and hence also for k, as the lower bound in (3.11) is 
only larger). However, this bound can be improved by choosing another 
and treating the case < 6 < 9o separately. To demonstrate this, take for 
instance Oq = ^. Since the above analysis gave that T > ^ whenever < |, 
it follows that 

r > / ^ , = — for all i < < I . 

- 160(1 -9o) 80 2 - - 4 

For < 60, we essentially repeat this analysis of T, only this time the 
respective value of 9 (that is, the maximum value it can attain) is ^. One 
can thus verify that in that case, 

(0t)3 

f{t, e) <et- ^-^ for 6* = i and any < t < 2.7 , 
6 

and the above argument then shows that 

t > — = — for ah < t < 2.7 . 
3 12 

On the other hand, 

9 11 
f > > — for alH > 2.7 , 

thus for 9 = ^ we have T > j2 for all t > 0. This converts into the lower 
bound r > ^ , thus completing the proof with a final value of k = . ■ 

Remark 3.10. Note that the only places where we used the fact that 9 < ^ 
are the lower bound on g{t, 9, Ci) in (3.11) and the analysis of T, as defined 
in (3.12). In both cases, we actually only need to have 9 < 9i for some 
constant ^1 < 1, whose precise value might affect the final value of k. 

Using the above lemma, we are now ready to obtain the final ingredient 
required for the proof of the recursion relation of incorporated in 

Lemma 3.11. This lemma provides a recursive bound on a quantity that 
resembles m^, where instead of integrating over xl, we integrate over f{xi). 

Lemma 3.11. Let f and D* he as in (3.4) and Lemma 3.3 respectively. 
There exists a universal constant k > so that for K = |(1 — 9)k we have 

/(4w.^(o-//(4Mft(o< „.(;;;,„j . 

Proof. Clearly, 

I f{xi)dQt{o- 1 /(4)^Q."(e) 

- f{xl))dQt{0 - I {f{xi) - f{xl))dQ-{0 . 
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Applying Lemma 3.7, we then obtain that 
f{xi)dQt{i)- j f{xi)dQ-{i) 

= CI I ifixi) - /«))(tanh ^ - tanh ^)dQ,{0 , 
and similarly, 

my = C* j {xl - x*)(^tanh y - tanh Y^dQ^(^) . 

Let 

F{x) = {fix) - fixD) (tanh(x/2) - tanh«/2)) , 
G{x) = (x - xl) (tanh(x/2) - tanh(x;/2)) , 

and define A to be the probability measure on M as: 

A(x) = ({e : x« = x}) . 

According to this definition, we have 

J F{xi)dQyiO = J F{x)dA , and J G{xi)dQy{C) = j G{x)dK , 
and thus, by the above arguments, 

f{4)dQUC) - [ f{4)dQy{C) = C: [ F{x)dA , 



m^ = C* j G{x)dA . (3.13) 

Furthermore, notice that by (3.7) and the fact that / is odd and increasing 
for X > 0, 



(tanh I - tanh 1^ 
and so 

Jf{xi)dQt{0- I f{4)dQy iO 

<2G:I / (^^) ( tanh ^ - tanh ^)dA. (3.14) 

In our next argument, we will estimate J G{x)dA and J G{x)dA according 
to the behavior of F and G about x*. Assume that x* > 0, and note that, 
although the case of x* < can be treated similarly, we claim that this 
assumption does not lose generality. Indeed, if x* < 0, one can consider the 
boundary condition of r' = — r, which would give the following by symmetry: 



X, 



xi , x;(-e) = -x«(o , QU-o = Qvio ■ 
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Therefore, as /(•) and tanh(-) are both odd functions, we have that / F(x)dA 
and / G{x)dA will not change under the modified boundary condition, and 
yet X*' > as required. 
Define 

/- = (-oo,<] , /+ = [x;,oo). 
First, consider the case where for either / = /"*' or / = /~ we have 

r J,F{x)dA>lJF{x)dA, 

\ JjG{x)dA> \jF{x)dA . ^ ' 

In this case, the following holds 

F{x)dA^ (^j G{x)dA^ <A(^j^F{x)dA^ (^j^G{x)dA 



<A j F{x)G{x)dA <A j F{x)G{x)dA 



where in the second line we applied the FKG-inequality, using the fact that 
both F and G are decreasing in I~ and increasing in . The last inequality 
followed from the fact that F and G are always non-negative. Note that 

J F{x)G{x)dA = I (/(x)-/(x:))(:E-x:)(tanh|-tanh^)'dA, 

and recall that Claim 3.9 applied onto tanh(x) (which indeed has an even 
non-negative derivative cosh~'^(x) that is decreasing in x > 0) gives 

tanh ^ — tanh ^ 2 tanh ^ — x > y . 

Noticing that each of the factors comprising F{x)G{x) has the same sign as 
that of (x — X*), and combining this with (3.7), it thus follows that 



F(x)dAj G{x)dA 

<16| /(^^)(x-x:)(tanh|-tanh^)tanh(^^)dA. 

(3.16) 

Second, consider the case where for and I~ as above, we have 
r i\^F{x)dA>\jF{x)dA, 

\jj.G{x)dA>lJG{x)dA. ^ ' 

The following definitions of F and G thus capture a significant contribution 
of F and G to J FdA and j GdA respectively: 

/ F{s)=F{xl + s) 

{ ^, , A ^, I , for any s > , (3.18) 
I G{s) = G{xl - s) 
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By further defining the probability measure A on [0, oo) to be 

A(s) = A« - s)l|,^o} + HK + s) for any s > , (3.19) 
we obtain that 

F{x)dA < 2 / F{x)dA < 2 / F{x)dA , 
Ji+ Jo 

/P /"OO 
G{x)dA <2 J G{x)dA <2 J G{x)dA . 

With both F and G being monotone increasing on [0,cxd), applying the 
FKG-inequality with respect to A now gives 

j F{x)dA^ (^j G{x)dA^ <A F{x)G{x)dA 
: 4 {fix: + s)- /«)) (^tanh - tanh 



i-s) (^tanh - tanh dA . 



Returning to the measure A, the last expression takes the form 
4/ (/(x)-/(x:))('tanh|-tanh^ 



• (x — xt) (tanh — — tanh — | dA 

^ "'^ V 2 2 J 

+4 if (2x1 -x)- fixD) (tanh - tanh ^) 

(X X* \ 

tanh — — tanh j dA . 



We now apply (3.7) and Claim 3.9 (while leaving the term (tanh | —tanh ^) 
unchanged in both integrals) to obtain that 



F{x)dA^ (^j G{x)dA 

<IQ j /(^^) (tanh I - tanh (x - <) tanh (^^^)dA . 

That is, we have obtained the same bound as in (3.16). 

It remains to deal with the third case where for and /~ as above, 

r j,.F{x)dA>\!F{x)dA, 

1 j,,G{x)dA>\!G{x)dA. ^"-"'^ 
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In this case, we modify the definition (3.18) of F and G appropriately: 



/ F{s) = F{xl + s) , , „ 

< _ for any s > , 

and let A remain the same, as given in (3.19). It then follows that 

/r /"OO 
F{x)dA <2 J F{x)dA <2 J F{x)dA , 

G{x)dA < 2 / G{x)dA < 2 / G{x)dA , 
Ji+ Jo 

with F and G monotone increasing on [0, oo). By the FKG-inequality, 

(^j F{x)dA^ (^j G{x)dA^ <A F{x)G{x)dA 

= Aj^ {fixl -s)- f{xl)) (tanh - tanh ^) 

• s (^tanh - tanh ^ j dA . (3.21) 

As before, we now switch back to A and infer from (3.7) and Claim 3.9 that 



F{x)dA^ (^j G{x)dA 
< 16 1 /(^^) (tanh I - tanh ^) (x - x^) tanh {^^)dA 



that is, (3.16) holds for each of the 3 possible cases (3.15), (3.17) and (3.20). 
Altogether, this implies that 



j f{xl)dQt{i) - I f{xi)dQ-{0 
= {C:f (I F{x)dA^ (1 G{x)dA^ 

<lQCf j /(^^) (tanh ^- tanh ^) _ tanh (^^)dA 
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Therefore, recalling (3.14) and choosing K = ^(1 — 9)k,, where n is as given 
in Lemma 3.8, we have 

/(4)dQ+(0- / f[xl)dQ~[i)]{l + Km,) 



/(^^) [l + 4^(1 - e)C:^^ tanh : 



• ( tanh tanh — ) dA 

V 2 2 ; 

< 2C: 1 (l/Z?:)^^^ ( tanh | - tanh ^)dA = Q^J G{x)dA . 

where the inequality in the last line is by Lemma 3.8 for 6 = \x — x*\/2 (the 
case X < X* follows once again from the fact that / is odd) and a choice of 
Ci = C* = 2cosh2(4/2) > 2 and C2 = [l/DD (recall that, by definition, 
1/D* = 1 + (^C* — 1)(1 — 0^) > 1, satisfying the requirements of the lemma). 
Therefore, (3.13) now implies that 

6m^ 



fi4)dQ:iO- f{xl)dQ-{i)< 



as required. 



Combining Lemmas 3.5 and 3.11, we deduce that there exists a universal 
constant k > such that 

m^< > , , , . (3.22 

The proof will now follow from a theorem of [31] , that links a function on 
the vertices of a tree T with its L^-capacity according to certain resistances. 

Theorem 3.12 ([31, Theorem 3.2] (reformulated)). Let T be a finite tree, 
and suppose that there exists some K > and positive constants {a^ : v G T} 
such that for every v €z T and x > 0, 

9v{x) < ayx/ (1 + Kx) . 

Then any solution to the system x^ = Ylw{v w)eT dwi^w) satisfies 

Xp < cap2(r)/K , 

where the resistances are given by R{u.v) = W{xy)&^{pv)^y'^ ' '^{p^'^) 
denoting the simple path between p and v. 

Together with inequality (3.22), the above theorem immediately gives 

cap2(f) 
K(l-e)/4 ' 

completing the proof of Theorem 3.2. ■ 
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Proof of Proposition 3.1. In order to obtain the required result from 
Theorem 3.2, recah the definition of xf, for v G T, according to which we 
can write 

fi^a{p) = l) = (l+ tanh(x«/2 + /i)) /2 , 

where h is the mentioned external field at the root p. By monotone coupling, 
we can construct a probability measure Qc on the space {(^, C') '■ C ^ C'} such 
that the two marginal distributions correspond to and Q~ respectively. 
It therefore follows that 

A = I {p^{a{p) = 1) - 0'{a{p) = l))dQ, 

= \ j (tanh(a;^/2 + h) - tanh(xf /2 + h)^dQc 

1 /■ 4 - 4 ^ ^ cap2(f ) 
< - / —dQc = -rrir, < — — , 

where the last inequality follows from Theorem 3.2 using the same value of 
^ ^ iQo- This completes the proof. ■ 



4. Upper bound on the inverse-gap and mixing time 

This section is devoted to the proof of the main theorem. Theorem 1, 
from which it follows that the mixing time of the continuous-time Glauber 
dynamics for the Ising model on a 5-ary tree (with any boundary condition) 
is poly-logarithmic in the tree size. 

Recalling the log-Sobolev results described in Section 2, it suffices to 
show an upper bound of 0(n log^'^^n) on inverse-gap of the discrete-time 
chain (equivalently, a lower bound on its gap), which would then imply an 
upper bound of 0{nlog^~^'^ n) for the mixing-time (and hence also for 
the total- variation mixing-time). 

The proof comprises several elements, and notably, uses a block dynam- 
ics in order to obtain the required upper bound inductively. Namely, we 
partition a tree on n vertices to blocks of size roughly n^~" each, for some 
small a > 0, and use an induction hypothesis that treats the worst case 
boundary condition. The main effort is then to establish a lower bound on 
the spectral-gap of the block dynamics (as opposed to each of its individual 
blocks). This is achieved by Theorem 4.1 (stated later), whose proof hinges 
on the spatial-mixing result of Section 3, combined with the Markov chain 
decomposition method. 

Throughout this section, let 6 > 2 be fixed, denote by Pc = arctanh(l/\/6) 
the critical inverse-temperature and let 6 = tanh /3c- 
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Figure 1. Block dynamics for the Ising model on the tree: 
illustration shows the distinguished block B(p, r) as well as 
a representative block of the form B{v,r) for v £ Hi. 

4.1. Block dynamics for the tree. In what follows, we describe our choice 
of blocks for the above mentioned block dynamics. Let h denote the height 
of our 6-ary tree (that is, there are b'^ leaves in the tree), and define 

£ = ah , r = h-e, (4.1) 

where < a < | is some (small) constant to be selected later. 

For any v £ T, let B(v, k) be the subtree of height k — 1 rooted at v, that 
is, B{v, k) consists of k levels (except when v is less than k levels away from 
the bottom of T). We further let denote the fc-th level of the tree T, 
that according to this notation contains b'^ vertices. 

Next, define the set of blocks B as: 

B = {B{v, r) : v£H(,U {p}} for £, r as above. (4.2) 

That is, each block is a 6-ary tree with r levels, where one of these blocks is 
rooted at p, and will be referred to as the distinguished block, whereas the 
others are rooted at the vertices of Hi. 

The following theorem establishes a lower bound on the spectral gap of 
the above-specified block dynamics (with blocks B). 

Theorem 4.1. Consider the Ising model on the b-ary tree at the critical 
inverse-temperature (3c and with an arbitrary boundary r. Let gap^ be the 
spectral gap of the corresponding block dynamics with blocks B as in (4.2). 
The following then holds: 

S^Pb > 4(6^V 1) " K{l-9){l-2a)) ' 
where k > is the absolute constant given in Theorem 3.2. 
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Given the above theorem, we can now derive a proof for the main result. 

4.2. Proof of Theorem 1. By definition, as 6 > 2, we have that 

^ = tanh/3.= i=<-L, 

hence we can readily choose an absolute constant < a < 1 such that 

^("^ = 8('-.(l-e)(l-2a)J>°- 

Let n/i = X]j=o ^ number of vertices in a 6-ary tree of height h 

excluding its leaves, and let gap^ be the spectral gap of the (single-site) 
discrete-time Glauber dynamics for the Ising model on a 6-ary tree of height 
h and boundary r (in the special case of a free boundary condition, rih 
should instead include the leaves). Further define 

gh = Uh min gap^; . 

T 

Recalling the definition of B according to the above choice of a, we have 
that each of its blocks is a tree of height r = (1 — a)/i, and that 

sup #{5 e ^ : X G B} = 2 , 

as each of the vertices in levels ^, ^ -|- 1, . . . , r is covered precisely twice in 
while every other vertex is covered precisely once. 

Hence, by Proposition 2.5 and Theorem 4.1, it now follows that for any 
/i > l/a (such that our choices of ^, r in (4.1) are both non-zero) we have 

Having established the induction step, we now observe that, as a is constant, 
clearly > c' holds for any k <l/a and some fixed c' = c'(a) > 0. Hence, 

g, > c'(c(a))^°^-'^(^/') = c'h-'<^)/'<^-^ , 

that is, there exists an absolute constant M (affected by our choice of the 
absolute constants k, a) so that the inverse-gap of the continuous-time dy- 
namics with an arbitrary boundary condition r is at most g^^ = 0{h^), as 
required. ■ 

4.3. Proof of Theorem 4.1. In order to obtain the desired lower bound 
on the spectral gap of the block dynamics, we will apply the method of 
decomposition of Markov chains, described in Subsection 2.6. To this end, 
we will partition our configuration according to the spins of the subset 

S = B{p,£-l) . 
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Note that S is strictly contained in the distinguished block B{p,r), and 
does not intersect any other B £ B. For rj G {±1}"^, denote the set of 
configurations which agree with rj by 

= {fj G O : 0-5 = r/} . 

Following the definitions in Subsection 2.6, we can now naturally decompose 
the block dynamics into a projection chain P on {±1}'^ and restriction chains 
Pri on $7^ for each rj G {±1}'^. With Theorem 2.6 in mind, we now need 
to provide suitable lower bounds on gap'^ and gap^, the respective spectral 
gaps of P and given the boundary condition r. 

We begin with the lower bound on the restriction chain gapj^, formulated 
in the next lemma. 

Lemma 4.2. For any boundary r and rj G {±1}'^, the spectral gap of the 
restriction chain on the space satisfies gapJJ > 1/(6^ + 1). 

Proof. Recall that the restriction chain P^ moves from cr G 17,; to c' G 17.^ 
(that is, o" and a' both agree with r] on S) according to the original law 
of the chain, and remains at a instead of moving to any a' ^ O^. By 
definition of our block dynamics, this means that with probability 6^/(6^ + 1) 
we apply a transition kernel Qi, that selects one of the blocks rooted at Hi 
to be updated according to its usual law (since S and all of these blocks are 
pairwise disjoint). On the other hand, with probability 1/(6^ + 1), we apply 
a transition kernel Q2 that updates the distinguished block, yet only allows 
updates that keep S unchanged (otherwise, the chain remains in place). 

We next claim that the update of the distinguished block can only in- 
crease the value of gap^. To see this, consider the chain P^, in which the 
distinguished block is never updated; that is, Q2 described above is replaced 
by the identity. Clearly, since each of the vertices of T \ 5 appears in (pre- 
cisely) one of the non-distinguished blocks, the stationary distribution of 
P^ is again ^r;r?) the Gibbs distribution with boundary conditions r] and r. 
Therefore, recalling the Dirichlet form (2.2), for any / we clearly have 

^p:M) = 2 ^ [fix)-fiy)]'f,^-''^{x)P;,{x,y) 

<l E [f{x)-f{y)ff^'''\x)Pr,{x,y)=£p^{f), 

and thus, by the spectral gap bound in terms of the Dirichlet form (2.1), 

gap(P^) > gap(P^) . (4.3) 

It remains to analyze the chain P^, which is in fact a product chain, and 
as such its eigenvalues can be directly expressed in terms of the eigenvalues 
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of its component chains. This well known fact is stated in the following 
straightforward claim (cf., e.g., [1, Chapter 4] and [22, Lemma 12.11]); we 
include its proof for completeness. 

Claim 4.3. For j £ [d], let Pj be a transition kernel on with eigenvalues 
Aj. Let I' be a probability distribution on [d], and define P' , the transition 
matrix of the product chain of the Pj-s on = fix x x • • • x by 

d 

P'{{xi,...,xd),{yi,...,yd)) = ^u{j)Pj{xj,yj) JJ 1 

j=l i-.ij^j 

Then P' has eigenvalues ^Yl'j=i '• Aj £ Aj| (with multiplicities) . 

Proof. Clearly, by induction it suffices to prove the lemma for d = 2. In this 
case, it is easy to verify that the transition kernel P can be written as 

P = j.(l)(Pi /f^J + z.(2)(/ni <S) P2) , 

where (8) denotes the matrix tensor- product. Thus, by tensor arithmetic, 
for any u,v, eigenvectors of Pi,P2 with corresponding eigenvalues Ai,A2 
respectively, {u<^v) is an eigenvector of P with a corresponding eigenvalue 
of z^(l)Ai + z^(2)A2, as required. ■ 

In our setting, first notice that Qi itself is a product chain, whose com- 
ponents are the b^ chains, uniformly selected, updating each of the non- 
distinguished blocks. By definition, a single block-update replaces the con- 
tents of the block with a sample according to the stationary distribution 
conditioned on its boundary. Therefore, each of the above mentioned com- 
ponent chains has a single eigenvalue of 1 whereas all its other eigenvalues 
are 0. 

It thus follows that (a lazy version of Qi) is another product chain, 
giving Qi probability 6^/(6^-1-1) and the identity chain probability l/(6^-|-l). 
By Claim 4.3, we conclude that the possible eigenvalues of P^ are precisely 

{iT^ + iT^i:;-^-^.^ {».!}} ■ 

In particular, gap(i-'^) = 1/(6^ -|- 1), and (4.3) now completes the proof. ■ 

It remains to provide a bound on gap"^, the spectral gap of the projection 
chain in the decomposition of the block dynamics according to S. This is 
the main part of our proof of the lower bound for the spectral gap of the 
block dynamics, on which the entire proof of Theorem 1 hinges. To obtain 
this bound, we relate the projection chain to the spatial-mixing properties 
of the critical Ising model on the tree under various boundary conditions, 
studied in Section 3. 
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Lemma 4.4. For any boundary t, the spectral gap of the projection chain 
P on the space {il}*^ satisfies 

^^P^ -¥T~i{^- .{1 - e)ii - 2a)) ' 

where n > is the absolute constant given in Proposition 3.1. 

We prove this lemma by establishing a certain contraction property of the 
projection chain P. Recall that -P(??, ?/'), for r/, jy' E {±1}'^, is the probability 
that completing i] into a state a according to the stationary distribution 
(with boundary r/ and r) and then applying the block dynamics transition, 
gives some a' that agrees with rj' on S. 

Let S* = Hi_i denote the bottom level of S, and notice that in the above 
definition of the transition kernel of P, the value of the spins in S\S* do not 
affect the transition probabilities. Therefore, the projection of the chain P 
onto S* is itself a Markov chain, which we denote by P* . In fact, we claim 
that the eigenvalues of P and those of P* are precisely the same (with the 
exception of additional 0-eigenvalues in P). To see this, first notice that the 
eigenfunctions of P* can be naturally extended into eigenfunctions of P with 
the same eigenvalues (as P* is a projection of P). Furthermore, whenever 
Vi V2 ^ S agree on S*, they have the same transition probabilities to any 
r]' £ S, thus contributing a 0-eigenvalue to P. It is then easy to see that all 
other eigenvalues of P (beyond those that originated from P*) must be 0. 
Altogether, 

gap(P*) = gap(P) ( = gap- ) , (4.4) 

and it remains to give a lower bound for gap(P*). The next lemma shows 
that P* is contracting with respect to Hamming distance on {±1}"^*. 

Lemma 4.5. Let and Y^* be instances of the chain P* , starting from ip 
and ijj respectively. Then there exists a coupling such that 

E,,dis.(xr,fn < (j^ + . s^i^i^) dist(,,v,) . 

Proof. Clearly, 'd = ip the lemma trivially holds via the identity coupling. 
In order to understand the setting when ip ^ ip, recall the definition of the 
chain P* , which has the following two possible types of moves Ei and E2 ■ 

(1) With probability 1 — the block dynamics updates one of the 

non-distinguished blocks: denote this event by Ei. Since this oper- 
ation does not affect the value of the spins in the subset S (and in 
particular, in S*), the projection chain P remains in place in this 
case (and so does P*). 
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(2) With probability j^rp[, the distinguished block is being updated: 
denote this event by £'2- By the discussion above, this is equivalent to 
the following. Let rj denote the current state of the chain P* . First, 
T\S is assigned values according to the stationary distribution with 
boundary rj and r. Then, the distinguished block B(p, r) is updated 
given all other spins in the tree, and the resulting value of S (and 
hence also of S*) is determined by the new state of the projection 
chain. 

By the triangle inequality, it suffices to consider the case of dist{(p,ip) = 1. 
Suppose therefore that (p and agree everywhere on S* except at some 
vertex g, and that without loss of generality, 

ifig) = 1 , V(£>) = -1 . 

Crucially, the above mentioned procedure for the event E2 is precisely 
captured by the spatial-mixing properties that were studied in Section 3. 
Namely, a spin of some site v G S* is propagated down the tree (with 
boundary condition r), and then the new value of S* is reconstructed from 
level r + 1, the external boundary of B(p,r). We construct a monotone 
coupling that will accomplish the required contraction property. 

First, when propagating the sites v £ S* with v ^ g, we use the identity 
coupling (recall that (p{v) = ip{v) for all v ^ g). Second, consider the 
process that the spin at g undergoes. For (p, a positive spin is propagated 
to Tg (with boundary condition r) and then reconstructed from level r + 1 
in the tree T (which corresponds to level r — i + 1 in the subtree Tg), with 
an additional boundary condition from T \Tg that translates into some 
external field. For ^/>, a negative spin is propagated analogously, and notice 
that in its reconstruction, the exact same external field applies (as T \ Tg 
was guaranteed to be the same for ip and ip). 

Therefore, applying Proposition 3.1 on the tree Tg with respect to the 
subtree T = B{g, r — I + 1), we can deduce that 

117 fv*f \ M F ^ ^ cap2(^(g,r -£ + 1)) 

^ip,i, [X^ [g) - [g) \E2j < _ , (4.5) 

where k > and the resistances are assigned as 

i?(„,,) = (tanh/3,)-2d-t(^''-) . 

We now turn to estimating the L'^-capacity, which is equivalent to the ef- 
fective conductance between g and dB{g, r — £ + 1). This will follow from 
the well-known Nash- Williams Criterion (cf., e.g., [24]). Here and in what 
follows, -Reff = l/C'eff denotes the effective resistance. 
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Lemma 4.6 (Nash- Williams Criterion [30]). // {11^}^^^ is a sequence of 
pairwise disjoint cutsets in a network G that separate a vertex v from some 
set A, then 

j eeUj 

In our case, G is the 6-ary tree -B(^), r — £ + 1), and it is natural to select 
its different levels as the cutsets Ilj. It then follows that 

i?eff(£'^5S(^,r-£+l)) > ^{h^Q'^^y^ = r-i^ \ , (4.6) 

fc=i 

where we used the fact that tanh /3c = ^ = 1/ \/^. It therefore follows that 

cap2(5(^,r -£+!)) < , 

which, together with (4.5), implies that 

E.,* - Ym I E,) < ^(j,;,,^,,) . (4.7) 

Unfortunately, asides from controlling the probability that the spin at g will 
coalesce in and ijj, we must also consider the probability that g would 
remain different, and that this difference might be propagated to other ver- 
tices in S* (as part of the update of B{p,r)). Assume therefore that the 
we updated the spin at g and indeed Xl(g) ^ Y-^{g), and next move on 
to updating the remaining vertices of S* . Since our propagation processes 
corresponding to X* and Y* gave every vertex in T \ the same spin, it 
follows that each vertex v £ S* , v ^ g, has the same external field in X* 
and Y* , with the exception of the effect of the spin at g. 

We may therefore apply the next lemma of [3], which guarantees that 
we can ignore this mentioned common external field when bounding the 
probability of propagating the difference in g. 

Lemma 4.7 ([3, Lemma 4.1]). Let T be a finite tree and let v ^ w be 
vertices in T. Let {Jg > : e G E{T)} be the interactions on T, and let 
{H{u) £ M : u € ^(r)} be an external field on the vertices of T. We 
consider the following conditional Gibbs measures: 

IJ^'^ : the Gibbs measure with external field H conditioned on cr{v) = 1. 
: the Gibbs measure with external field H conditioned on cr{v) = —1. 

Then iJ.'^'^ {a{w)) — fi~'^ {a{w)) achieves its maximum at H = 0. 
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In light of the discussion above, Lemma 4.7 gives that 



IY.{XI{v)-Y*{v))\E2 



ves* 



k=l 

^ l + {b-l){i-l)/b ^ l + {b-l)£ 
- K{l-e){r-i) ~ bK{l - e){r - I) ' 

where in the first inequahty we used the propagation property of the Ising 
model on the tree (Claim 2.1), and in the second one we used the fact that 
6 = tanh(/9c) = as well as the estimate in (4.7). 

We conclude that there exists a monotone coupling of and with 

l + (6-l)^ 



^-H^-*(^-^^*)I^^J^Mi-^)(.-.)' 

which then directly gives that 

1 l + (6-l)£ 



E^,^(dist(Xi*,yi*)) < 7T-T + 



+ 1 ¥ + 1 bK{l-e){r-l) ' 
as required. ■ 

The above contraction property will now readily infer the required bound 
for the spectral gap of P* (and hence also for gap"^). 

Proof of Lemma 4.4. The following lemma of Chen [8] relates the con- 
traction of the chain with its spectral gap: 

Lemma 4.8 ([8]). Let P be a transition kernel for a Markov chain on a 
metric space il. Suppose there exists a constant l such that for each x,y £ 0,, 
there is a coupling (Xi,Yi) of P(x,-) and P{y,-) satisfying 

E^,j,(dist(Xi,yi)) < i dist(x,y) . (4.8) 

Then the spectral gap of P satisfies gap > 1 — l. 

By Lemma 4.5, the requirement (4.8) is satisfied with 

_ b^ 1 l + (b-l)^ 

~ ¥TT ^ fe^TT " 6k(i - e){r - 1) ' 

and hence 



+ lV bK{l-e){r-£) 
- ¥+l{^ ~ At(l-e)(l-2a),' ' ^^"^^ 
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where in the last inequahty we increased 1 + (6 — 1)£ into bi to simphfy the 
final expression. This lower bound on gap(P*) translates via (4.4) into the 
desired lower bound on the spectral gap of the projection chain, gap"^. ■ 

We are now ready to provide a lower bound on the spectral gap of the 
block dynamics, gap^, and thereby conclude the proof of Theorem 4.1. By 
applying Theorem 2.6 to our decomposition of the block dynamics chain Pg, 

„-=r,^ ^ sap gap • gap^^^i^ 

gapg > — A — — - — z — , (4.10) 
" 3 37 + gap 

where 

gaPmin= min,gap^> 7= max niax ^ Pb^x^v) ■ 

Lemma 4.2 gives that gapj^j^^ > 1/(6^ + 1), and clearly, as the spins in S can 
only change if the distinguished block is updated, 7 < 1 /(6^ + 1). Combining 
these two inequalities, we obtain that 

gap • gap^,, _ gap_ ^ / _ > \ f ^ ^ gap) (4.11) 



37 + gap 1 + 37/gap - (6^ + 1) + 3/gap " 4 Vft^ + 1 
with room to spare. Together with (4.10), this implies that 

g-PS > A ^gap , 

and Lemma 4.4 now gives that 

gapg > .... , 1 



+ 1)V k(1 -6')(1 -2a) 

as required. This concludes the proof of Theorem 4.1, and completes the 
proof of the upper bound on the mixing time. ■ 

Remark 4.9. Throughout the proof of Theorem 1 we modified some of the 
constants (e.g., (4.9), (4.11) etc.) in order to simplify the final expressions 
obtained. By doing the calculations (slightly) more carefully, one can obtain 
an absolute constant of about 300 for the upper bound in Theorem 1. 



5. Lower bounds on the mixing time and inverse-gap 

In this section, we prove Theorem 3, which provides lower bounds on 
the inverse-gap and mixing time of the critical Ising model on a 6-ary tree 
with free boundary. Throughout this section, let 6 > 2 be fixed, and set 
e = tanh^c = 4r- 
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5.1. Lower bound on the inverse-gap. The reqmred lower bound will 
be obtained by an application of the Dirichlet form (2.1), using a certain 
weighted sum of the spins as the corresponding test function. 

Proof of Theorem 3, inequality (1.2). Let T be a 6-ary tree, rooted at 
p, with h levels (and n = X^^^g vertices). We will show that 

gap > -^rih . 

For simplicity, we use the abbreviation d{v) = dist(p, f), and define 
g{a) = ^9<'"^a{v) for a G 17 . 

By the Dirichlet form (2.1), and since P{cr, a') < ^ for any a,a' £ in the 
discrete-time dynamics, we have that 



<T,a' 

< \ u.^^Y.^g{a) - g{a')f^,{<y)P{a,a') < ^ -(2e'=)2 < 



- Afc 2h 



2 o- ^-^ 2 ^-^ n n 

a' k=0 

On the other hand, the variance of g can be estimated as follows. 

Var^5 = Var^ (^0'^(^)a(^;)) = ^ Cov^(a(n), a(u')) 

= E ^'^"^+'^"'^C0V^(^(t^),^(^))l{.A«,=.} ' 
u,v,w£T 

where the notation (n A w) denotes their most immediate common ancestor 
(i.e., their common ancestor z with the largest d{z)). Notice that for each 
V £ T, the number of u,w that are of distance i,j from v respectively and 
have V = u Aw IS precisely U ■ {h — 1)6^~^, since determining u immediately 
rules 6^"^ candidates for w. Furthermore, by Claim 2.1 we have 

Cow ^{a{u), a iw)) = ^'^W+'^W-s^'M ^ 

and so 

h h—kh—k 

= Yb^YY^ b"{b - l)^^?■-l^2fc+^+j^^+i 

k=0 i=0 j=0 

6-lA.. .,o 5-1 6-1 



k=0 
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Figure 2. Speed-up dynamics for the Ising model on the tree. 



Altogether, we can conclude that 

8{g) _ 6b 1 



gap < 



Var^ g 6—1 n/i^ ' 



(5.1) 



as required. 



5.2. Lower bound on the mixing-time. In order to obtain the required 
lower bound on the mixing time, we consider a "speed-up" version of the 
dynamics, namely a custom block-dynamics comprising a mixture of single- 
tons and large subtrees. We will show that, even for this faster version of 
the dynamics, the mixing time has order at least n log^ n. 

Let T be a 6-ary tree with h levels (and n = ^^=0 vertices). Consider 
two integers 1 < £ < r < /i, to be specified later. For every v G H^, 
select one of its descendants in Hj- arbitrarily, and denote it by w^. Write 
W = {wy : V £ Hi} as the set of all such vertices. Further define 

By = {T-u \ TwJ U {wv} (for each v G He) . 

The speed-up dynamics, (Xt), is precisely the block-dynamics with respect 
to 

B = {By:veHi}U U {u} . 

In other words, the transition rule of the speed-up dynamics is the following: 

(i) Select a vertex u £ V(T) uniformly at random. 

(ii) li u ^ W, update this site according to the usual rule of the Glauber 
dynamics. 

(iii) Otherwise, update given the rest of the spins, where v £ He is the 
unique vertex with u = w^. 
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Figure 3. Speed-up dynamics on the forest F and the sub- 
forest G. 



The fohowing theorem of [33] guarantees that, starting from ah-plus con- 
figuration, the speed-up Glauber dynamics indeed mixes faster than the 
original one. In what follows, write ^ v ii fi stochastically dominates ly. 

Theorem 5.1 ([33] and also see [32, Theorem 16.5]). Let ($7,5, vr) be 
a monotone system and let ^ be the distribution on 0, which results from 
successive updates at sites vi, . . . ,Vm, beginning at the top configuration. 
Define v similarly but with updates only at a subsequence Vi^, . . . ,Vi^. Then 
li^v, and ll/i — vtIItv ^ ||2^~'''"||tv- Moreover, this also holds if the sequence 
vi, . . . ,Vm and the subsequence ii, . . . , are chosen at random according to 
any prescribed distribution. 

To see that indeed the speed-up dynamics Xt is at least as fast as the 
usual dynamics, first note that any vertex u ^ W is updated according to 
the original rule of the Glauber dynamics. Second, instead of updating the 
block B^, we can simulate this operation by initially updating Wv (given its 
neighbors), and then performing sufficiently many single-site updates in B^. 
This approximates the speed-up dynamics arbitrarily well, and comprises a 
superset of the single-site updates of the usual dynamics. The above theorem 
thus completes this argument. 

It remains to estimate the mixing time of the speed-up dynamics Xf. To 
this end, define another set of blocks as follows: for every v £ Hi, let Ly 
denote the simple path between v and (inclusive), define the forest 

and put 

BF = {L,:ve Hi} U \J {u} . 

ueF\w 

We define Yt, the speed-up dynamics on F, to be the block-dynamics with 
respect to Bp above. This should not be confused with running a dynamics 
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on a subset of T with a boundary condition of the remaining vertices; rather 
than that, Yt should be thought of as a dynamics on a separate graph F, 
wliich is endowed with a natural onc-to-onc mapping to the vertices of T. 
Further note that, except for the singleton blocks in B, every block By E B 
in the block-dynamics Xf has a counterpart Ly C By 'm Yf. 

The next lemma compares the continuous-time versions of Xf and Yf 
(where each block is updated at rate 1), and shows that on a certain subset 
of the vertices, they typically remain the same for a substantial amount of 
time. 

Lemma 5.2. Let (Xf) and (!<) be the continuous-time speed-up dynamics 
on T and F respectively, as defined above. Let G = \^y^H^ Tyj^ and define 

T = \ui{Xt{u) 7^ Yt{u) for some u G V{G)} . 

Then there exists a coupling of Xt and Yt such that 

P(t >t)> expi-e'-^bh) . 

Proof. For two configurations a G {±1}"^ and r) G {±1}^, denote their 
Hamming distance on F by 

dist((7,77) = ^{a{v)j^riiv)} ■ 

The coupling of Xt and Yt up to time r can be constructed as follows: 

(1) Whenever a singleton block {u} with u eT\F is being updated in 

Xt, the chain Yt remains in place. 

(2) Otherwise, when a block B is updated in Xt, we update B f] F (the 
unique B' G Bp with B' C B) in Yf so as to minimize dist(Xt, Yf). 

For any w G W, define the stopping time 

= inf{t : Xt{w) ^ Yt{w)} , 

and notice that in the above defined coupling we have r = min„,gvF Tw, since 
W separates G\W from F. 

Let V G Hi and w = w,, G W, and suppose that block By is to be updated 
at time t < in Xt, and hence, as defined above, Ly is to be updated in 
Yt. By definition, at this time these two blocks have the same boundary 
except for at v, where there is a boundary condition in T (the parent of v) 
and none in F (recall v is the root of one of the trees in F) . 

We now wish to give an upper bound on the probability that this update 
will result in Xt{w) / Yt{w). By the monotonicity of the Ising model, it 
suffices to give an upper bound for this event in the case where v has some 
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parent z in F, and Xti^z) ^ yt{z). In this case, we can bound the probabihty 
that Xt{w) 7^ Yt{w) (in the maximal couphng) by an expression of the form 

^(/x+'^(aH)-/.-'^(aH)) 

as described in Lemma 4.7, where the external field H corresponds to the 
value of the spins in Tyj\{w}. Lemma 4.7 then allows us to omit the external 
field H at w, translating the problem into estimating the probability that a 
difference propagates from v to w. By Claim 2.1, we deduce that 

F{Xt{w)^Yt{w))<e'-' , 



and therefore 



Using the fact \W\ = b^, it follows that 



< r) = P t < min > exp{-0''-Vt) 



as required. 



With the above estimate on the probability that Xt and Yt are equal on 
the subgraph G up to a certain time-point, we can now proceed to studying 
the projection of Xt on G via that of Yt (being a product chain, Yj is much 
simpler to analyze). 

To be precise, let Xt and Yt denote the respective projections of Xt and 
Yt onto G, which as a reminder is the union of all trees T^„. Notice that 
Yt is precisely the continuous-time single-site Glauber dynamics on G, since 
the block update of in F translates simply into the single-site update of 
in G. On the other hand, Xt is not even necessarily a Markov chain. We 
next prove a lower bound on the mixing time of the Markov chain Yt. 

Lemma 5.3. Let Ht be the transition kernel of Yt, and let hg denote its 
corresponding stationary measure. Let gap' denote the spectral-gap of the 
continuous-time single-site dynamics on a b-ary tree of height h — r. Then 

\\Ht{l, •) - /icIlTV > I for any t < , 
where 1 denotes the all-plus configuration. 

Proof. Let T' denote a 6-ary tree of height h — r and n' vertices. Let P' be 
the transition kernel of the corresponding discrete-time single-site Glauber 
dynamics on T', let Hj^ be the transition kernel of the continuous-time version 
of this dynamics, and let fj,' be their corresponding stationary measure. 

By definition of G as a disjoint union of b^ copies of T', clearly is a 
product of b^ copies of identical and independent component chains on T'. 
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We can therefore reduce the analysis of Yt into that of Hj-, where the second 
eigenvalue of of its discrete-time counterpart P' plays a useful role. 

The following lemma ensures that P' has an increasing eigenfunction cor- 
responding to its second largest eigenvalue A'. 

Lemma 5.4 ([29, Lemma 3]). The second eigenvalue of the discrete-time 
Glauber dynamics for the Ising model has an increasing eigenfunction. 

Since the eigenspace of A' has an increasing eigenfunction, it also con- 
tains a monotone eigenfunction / such that |/(1)| = ||/||oo- Therefore, the 
transition kernel of the continuous-time chain satisfies 



H'tf) a) = (E«"*"'^(^')V)(i) 



e~^" 



E ^^/(i) = e-"'(^-^')V(l) . (5.2) 



k=0 

Since J fdfi' = 0, we have that 

mf){l)\ = I {H^a, y)f{y) - f{y)fi'{y)) \ < 2||/||oo||i^*(l, •) - m'IItv • 
y 

Plugging in (5.2) and using the fact that |/(1)| = ||/||oo, it follows that 

||i^;(l,-)-^'||TV>^e-"'(i-^')*. (5.3) 

In order to relate the product chain Yj to its component chain Y^, we will 
consider the Hellinger distance between certain distributions, defined next 
(for further details, see, e.g., [21]). First, define the Hellinger integral (also 
known as the Hellinger affinity) of two distribution /x and i/ on Q to be 

InifJ; '^) = Y1 V/^(x)zy(x) . 

The Hellinger distance is now defined as 

dnifJ', = \/2 - 2I-h(//,z/) . 
Clearly, for any two distributions fj, and v, 

and so d-^-i provides the following lower bound on the total variation distance: 

IIm - i'IItv > 1 - In{l^, ^) = ^dliifJ', v) ■ (5.4) 

Furthermore, the Hellinger distance also provides an upper bound on d^y, 
as the next simple inequality (e.g., [15, Lemma 4.2 (i)]) shows: 

IIm - '^Wtv < dnifJ-,!^) ■ (5.5) 
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To justify this choice of a distance when working with product chains, notice 
that any two product measures ^ = 

nr=i A*^'^ and V = nr=i satisfy 

n 

Iniy^.v) =Wln{i^^^ M"^) . (5.6) 
1=1 

Next, we consider the Hehinger integral of our component chains Indeed, 
combining the definition of d-n with (5.5), we get that 

/„(i?;(i,.),/i') < 1 - ^ii^;(i,-) -/x'lilv < 1 - ^e-2"'(i-^')* , 

where the last inequality is by (5.3). Therefore, applying (5.6) to the product 
chain B.t (the product of copies of i/^), we can now deduce that 

/7.(^*(l,-),^G)<(^l-V2(l-^')*n''' 
At this point, (5.4) gives that 

-2{\-\')tn' 

ll^t(i,-)-/"G||TV>l-(l ^ ) . 

Recall that by definition, gap' is the spectral-gap of the continuous-time 
version of P', and so gap' = n'(l — A'). Hence, if 

^log6-2 
- 2gap' 

then 

ll^td, •) - /^gIItv > 1 - exp ( - eVs) > ^ , 
as required. ■ 

The final ingredient required is the comparison between (the Gibbs 
distribution on G), and the projection of \i (the Gibbs distribution for T) 
onto the graph G. The following lemma provides an upper bound on the 
total-variation distance between these two measures. 

Lemma 5.5. Let [i and fic be the Gibbs distributions for T and G resp., 
and let jl denote the projection of ^ onto G, that is: 

m = ^({^ G {±1}"^ : (TG = ri}) (forrje {±1}« ) . 
Then ll/XG -/uIItv < b'^^9^^^-^\ 

Proof. Recalling that G is a disjoint union of trees {Ty^ : w G W}, clearly the 
configurations of these trees are independent according to fiQ- On the other 
hand, with respect to fl, these configurations are correlated through their 
first (bottom-most) common ancestor. Further notice that, by definition, 
the distance between Wi wj £ W in T is at least 2(r — i + 1), as they 
belong to subtrees of distinct vertices in H^. 
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To bound the effect of the above mentioned correlation, we construct a 
coupHng between fic and jl iteratively on the trees {T^ : w £ W}, generating 
the corresponding configurations t] and ij, as follows. Order W arbitrarily as 
W = {wi, . . . , Wfji}, and begin by coupling hg and jl on T^^ via the identity 
coupling. Now, given a coupling on Uj<fcr^j, we extend the coupling to 
r^j, using a maximal coupling. Indeed, by essentially the same reasoning 
used for the coupling of the processes Xt and 1^ on G in Lemma 5.2, the 
probability that some already determined Wi (for i < k) would affect Wk is 
at most Summing these probabilities, we have that 

IP {VT^, + ^T.J =^{j]{wu) / ^{Wk)) <{k- 

Altogether, taking another union bound over all k G [6^], we conclude that 

||/iG-/i||TV<IP(??/r?)<6V'-^) , 

completing the proof. ■ 

We are now ready to prove the required lower bound on tMix- 

Proof of Theorem 3, inequality (1.3). As we have argued above (see 
Theorem 5.1 and the explanation thereafter), it suffices to establish a lower 
bound on the mixing time of the speed-up dynamics Xt on T. By considering 
the projection of this chain onto G, we have that 

\\^i{Xt e •) - AillTv > l|IPi(^i G •) - AIItv , 
and recalling the definition of r as mit{{Xt)G 7^ iYt)G}i 

\\F,{Xt G •) - Mtv > W^iiYt G •) - AIItv - IP(t < t) 

> \\Fi{Yt G •) - /^gIItv - Hr <t)- WfiG - AIItv • 

Let gap and gap' denote the spectral-gaps of the continuous-time single-site 
dynamics on a b-aiy tree with h levels and h — r levels respectively (and free 
boundary condition), and choose t such that 

t < ""^"T^ . (5.7) 

Applying Lemmas 5.2, 5.3 and 5.5, we obtain that 

\\¥i{Xt G •) - MIItV > ^ - (l - expi-e'^-Vt)'^ - 62^e2(r-£) _ (5_g^ 

Now, selecting 
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and recalling that bO^ = 1, we have that the last two terms in (5.8) both 
tend to as /i — > oo, and so 

e-)-^l|Tv> • 

5 

In particular, for a sufficiently large h, this distance is at least 1/e, hence by 
definition the continuous-time dynamics satisfies tMix ^ t. We can can now 
plug in our estimates for gap' to obtain the required lower bounds on tMix* 
First, recall that by (5.1), 

66 1 

and so the following choice of t satisfies (5.7): 

It follows that the mixing-time of the continuous-time dynamics satisfies 

and the natural translation of this lower bound into the discrete-time version 
of the dynamics yields the lower bound in (1.3). 

Second, let g{h) be the continuous-time inverse-gap of the dynamics on 
the 6-ary tree of height h with free-boundary condition, and recall that by 
Theorem 1, we have that g is polynomial in h. In particular, 

g{h) < Cg{h/5) for some fixed C > and all h. 

Since by definition (gap')"-*^ = g{h — r) = g{h/b) and gap^-*^ = g{h)^ we can 
choose t to be the right-hand-side of (5.7) and obtain that for any large h 

^Mix > t > C'gap~^/i for some C" > fixed. 

Clearly, this statement also holds when both tuix and gap correspond to the 
discrete-time version of the dynamics, completing the proof. ■ 

6. Phase transition to polynomial mixing 

This section contains the proof of Theorem 2, which addresses the near 
critical Ising model on the tree, and namely, the transition of its (continuous- 
time) inverse-gap and mixing-time from polynomial to exponential in the 
tree-height. Theorem 2 will follow directly from the next theorem: 

Theorem 6.1. Fix b > 2, let e = e{h) satisfy < e < Eq for a suitably 
small constant eq, and let j3 = arctanh (y^ (1 -|- e) /6) . The following holds 
for the continuous-time Glauber dynamics for the Ising model on the b-ary 
tree with h levels at the inverse-temperature (3: 
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(i) For some ci > fixed, the dynamics with free boundary satisfies 

gap-i > ci ((1/e) A hf{l + e)'' . (6.1) 
(a) For some absolute constant C2 > and any boundary condition t 

gap-i <tMix<e^^(^''+^°s'') . (6.2) 

Throughout this section, let 6 > 2 be some fixed integer, and let T be 
a 6-ary tree with height h and n vertices. Define 9 = y^(1 + e) /b, where 
e = e(n) satisfies < e < eo (for some suitably small constant Eq < g to be 
later specified), and as usual write (3 = arctanh(0). 

Proof of Theorem 6.1. The proof follows the same arguments of the proof 
of Theorems I and 3. Namely, the upper bound uses an inductive step using 
a similar block dynamics, and the decomposition of this chain to establish 
a bound on its gap (as in Section 4) via the spatial mixing properties of the 
Ising model on the tree (studied in Section 3). The lower bound will again 
follow from the Dirichlet form, using a testing function analogous to the one 
used in Section 5. As most of the arguments carry to the new regime of /3 
in a straightforward manner, we will only specify the main adjustments one 
needs to make in order to extend Theorems I and 3 to obtain Theorem 6.1. 

Upper bound on the inverse-gap. Let < k < I be the universal 
constant that was introduced in Lemma 3.8 (and appears in Proposition 3.1 
and Theorem 3.2), and define 

A K 1 

"° - 20 - 8 • 

As 6 > 2 and e < eo ^ g) we have that ^ < |, hence Proposition 3.1 and 
Theorem 3.2 both hold in this supercritical setting. It therefore remains to 
extend the arguments in Section 4 (that use Proposition 3.1 as one of the 
ingredients in the proof of the upper bound on gap~^) to this new regime 
of p. 

Begin by defining the same block dynamics as in (4.2), only with respect 
to the following choice of £ and r (replacing their definition (4.1)): 

a = eo = k/20 , (6.3) 

£ = a[{l/e) Ah] , r = h-i . (6.4) 

Following the same notations of Section 4, we now need to revisit the argu- 
ments of Lemma 4.5, and extend them to the new value of = tanh /? = 
V(l + e)/fe. This comprises the following two elements: 

(1) Bounding the L^-capacity ca'P2{B{g,r — i)). 

(2) Estimating the probability that a difference in one spin would prop- 
agate to other spins, when coupling two instances of the chain P* . 
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Recalling the Nash- Williams Criterion (Lemma 4.6) and its application 
in inequality (4.6), the effective resistance between q and dB{g,r — i) is at 
least 

r-i+l r-£+l 

J2 {b'e^''y'= ^ (1+ £)-'= = i(i-(i + e)-('-^+i)) , 

k=l k=l 

which implies that 

cap2 {B{e, r -£ + !))< ^ _ ^^^^_(,_,) ■ (6.5) 

Now, if e > we have 

1 - (1 + e)-('-^) = 1 - (1 + ^)~ih-'2a/e) > ^ _ + ^yil~2a)/e 

> 1 - e-(i-2") > , 

- 2 ' 

2 

where the last inequality uses the fact that exp(— x) < 1 — x + ^ and that 
a > 0. Similarly, if e < then 

1 _ (1 + = 1 - (1 + e)-Ml-2«) > 1 _ ^-ehil-2a) 

, (eh(l-2a)f I - 2a 

> eh{l - 2a) - ^ > eh^- , 

where in the last inequality we plugged in the fact that eh < 1. Combining 
the last two equations with (6.5), we deduce that 

1 — 2a 

Using (4.5), it then follows that 

2{e V (l/h)) 



E^,^[Xlig)-Y{{e)\E2) < 



K(l-6')(l-2a) ■ 

By repeating the next arguments of Lemma 4.5 (without any additional 
essential changes), we obtain that under the monotone coupling, 

2(e V {l/h)) fl (1 + e)"[(i/£)Ah] _ i 



k{1 - e){l - 2a)\b 



^ 2(e V {l/h)) (l + e)°[(V^)^'»] - 1 ^ 2(e V {l/h)) e^I^^^'^l - 1 



< 



K{l-e){l-2a) e - K{l-e){l-2a) 

2(e V {l/h)) 2a[l A eh] 4a 



K{l-9){l-2a) e k{1 - e){l - 2a) ' 

where in the last line we used the fact that e^ — 1 < 2x for all < x < 1. 
Again defining gf^ = rih min,- gapj^, we note that all the remaining arguments 
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in Section 4 apply in our case without requiring any modifications, hence 
the following recursion holds for gh'. 

9h > c{a)gr = c(a)c/,,_a[(i/£)^/j] , (6.6) 

where 

^("^ = ^('-.(i-^Ki-2«)) • 

Recalling the definition (6.3) of a, since ^ < | and k < 1 we have that 

4a 2 8 

< 



K{l-e){l-2a) (1-0)(1O-k) 9' 
and so c(a) > 0. We now apply the next recursion over g^^.'- 

, ^ ^ fhk-{a/e) iihk>{l/e), 

[{l-a)hk ifhk<{l/e). 

Notice that by our definition (6.3), we have e < Eq = a. With this in mind, 
definition (6.4) now implies that for any h > 1/a we have £,r > 1. Thus, 
letting K = min{A; : hk < we can conclude from (6.6) that 

9hk > c{a)gh^^^ for all k < K, and hence 

gh > {c{a))^ ghji ■ 
By the definitions of and clearly 

K<^h + logi/(i_„) {h A (l/e)) = 0{eh + log h) . 

Since < clearly g^j^. > c' for some constant c' = c'{a) > 0, giving 

5h >c'(c(a))^ >e-^'^(^'^+'°s^) 

for some constant M = M(a) > and any sufficiently large n. By definition 
of gh, this provides an upper bound on gap^-*^, and as tMix = O (gap^-*^ log^ n) 
(see Corollary 2.4 in Section 2), we obtain the upper bound on ^mix that 
appears in (6.2). 

Lower bound on the inverse-gap. We now turn to establishing a lower 
bound on the inverse-gap. Define the test function g to be the same one 
given in Subsection 5.1: 

g(a) = ^0'ii^*(^'^)a(^;) . 

By the same calculations as in the proof of Theorem 3 (Subsection 5.1), we 
have that 

£ig) < \ X: = - E(l + = 2 (l±£^!i^ , (6.7) 

2 ^-^ n n ^-^ n e 

k=0 k=0 
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whereas 

Var,(,) = V E ( E ^^^") ^ = V Ed + Ed + 

A:=0 j=0 fc=0 j=0 



, ^ h /i— fe „ , ^ h h—k 

- 1 ^ ,l..9i,/^ ,,„oA2 6-1 



5-1 



A:=0 
ft 



^ ((1 + e)2^-^+2 _ 2(1 + e)'^+i + (1 + e)/^) 



A:=0 

= ^ ((1 + 6)2'^+' - (2/i + 3) e(l + e)'^+i - l) . (6.8) 
When e > 8/h we have 

1(1 + _ (2/, + 3)e > |(/, + 2) + ^ + ^) - (2/i + 3)e 

>(/. + 2)eQ+e^-2) >4, 
and therefore in this case (6.8) gives 

Var,(<;) > -^^-^l . (6.9) 

Combining (6.7) and (6.9), the Dirichlet form (2.1) now gives that 

4b £2 

6 - ln(l + £ 

On the other hand, when < e < 8/h we stih have bO^ > 1 and hence 



gap < ^ ^ , fore>8//i. (6.10) 



Var,(,) = '-^i:b'e^'(Xb'e^'y > VE(^-^)' ^ 

k=0 i=0 k=0 

In addition, using the fact that the expression [(1 + e)^+^ — l]/e in (6.7) is 
monotone increasing in e, in this case we have 

, 2 (l + (8//i))^+^ -1 7, , 

where the last inequality holds for any h > 20. Altogether, the Dirichlet 
form (2.1) yields (for such values of h) 

3e^6 1 

g^P^53i;^ ^or0<e<8/h. (6.11) 
Combining (6.10) and (6.11), we conclude that 



gap 



< 



Se^^b r . .1-1 



5-1 



n(l + e)'^((l/e) A hy 
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where we used the fact that (1 + e)^ < e. This gives the lower bound on 
gap^^ that appears in (6.1), completing the proof of Theorem 6.1. ■ 

7. Concluding remarks and open problems 

• We have established that in the continuous-time Glauber dynamics for 
the critical Ising model on a regular tree with arbitrary boundary con- 
dition, both the inverse-gap and the mixing-time are polynomial in the 
tree-height h. This completes the picture for the phase-transition of the 
inverse-gap (bounded at high temperatures, polynomial at criticality 
and exponential at low temperatures), as conjectured by the physicists 
for lattices. Moreover, this provides the first proof of this phenomenon 
for any underlying geometry other than the complete graph. 

• In addition, we studied the near-critical behavior of the inverse-gap 
and mixing-time. Our results yield the critical exponent of /3 — (3c, as 
well as pinpoint the threshold at which these parameters cease to be 
polynomial in the height. 

• For further study, it would now be interesting to determine the precise 
power of h in the order of each the parameters gap^^ and tmix at the 
critical temperature. In the free-boundary case, our lower bounds for 
these parameters in Theorem 3 provide candidates for these exponents: 

Question 7.1. Fix b > 2 and let /Sc = arctanh(l/\/6) be the critical 
inverse-temperature for the Ising model on a b-ary tree of height h. Does 
the corresponding continuous-time Glauber dynamics with free boundary 
condition satisfy gap^-*^ x and ^mix x /i^ ? 

• Both at critical and at near-critical temperatures, our upper bounds for 
the inverse-gap and mixing-time under an arbitrary boundary condition 
matched the behavior in the free-boundary case. This suggests that 
a boundary condition can only accelerate the mixing of the dynamics, 
and is further supported by the behavior of the model under the all- plus 
boundary, as established in [27]. We therefore conjecture the following 
monotonicity of gap~^ and t^ix with respect to the boundary condition: 

Conjecture 7.2. Fix b >2 and P > 0, and consider the Ising model on 
a b-ary tree with parameter (3. Denote by gap and fmix the spectral-gap 
and mixing time for the Glauber dynamics with free boundary, and by 
gap"^ and t1^^^ those with boundary condition r. Then 

gap < gap'^ and tuix > ^mix for any r. 

• A related statement was proved in [26] for two-dimensional lattices at 
low temperature: It was shown that, in that setting, the spectral-gap 
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under the all-plus boundary condition is substantially larger than the 
spectral-gap under the free boundary condition. In light of this, it would 
be interesting to verify whether the monotonicity property, described in 
Conjecture 7.2, holds for the Ising model on an arbitrary finite graph. 
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